LEE- YANG PROBLEMS AND THE GEOMETRY OF 
MULTIVARIATE POLYNOMIALS 



JULIUS BORCEA AND PETTER BRANDEN 

Abstract. We describe all linear operators on spaces of multivariate polyno- 
mials preserving the property of being non- vanishing in open circular domains. 
This completes the multivariate generalization of the classification program 
initiated by Polya-Schur for univariate real polynomials and provides a natu- 
ral framework for dealing in a uniform way with Lee- Yang type problems in 
statistical mechanics, combinatorics, and geometric function theory. 



1. Introduction 

In O I12j Lee and Yang proposed the program of analyzing phase transitions in 
terms of zeros of partition functions and proved a celebrated theorem that may be 
stated as follows. Recall that the partition function of the Ising model (at inverse 
temperature 1) may be written as 

z(h 1 ,...,h n )= 

o-e{-i,i}» 

where a ■ h = Y^i=i / i ( fJ ) = e^-' i -'= 1 J ^ <7t ' T ^ ; the Jy are coupling constants and 
the hi are external (magnetic) fields sometimes also called fugacities. 

Theorem 1.1 (Lee- Yang [6]). If Jy > for all 1 < i,j < n then 

(a) Z(li\, . . . , h n ) ^ whenever Re(/ij) > 0, 1 < i < n; 

(b) All zeros of Z{h, . . . , h) lie on the imaginary axis. 

The important consequence of Theorem ll.il (b) is that the zeros of the partition 
function of the Ising model in the ferromagnetic regime (i.e., when all coupling con- 
stants are non-negative) accumulate on the imaginary axis in the complex fugacity 
plane and a phase transition occurs only at zero magnetic field. 

In their proof of Theorem 1 1.1 1 Lee and Yang used a theorem of Polya on the zero 
distribution of entire functions, thus establishing a connection between statistical 
mechanics and a classical topic in geometric function theory. This connection has 
become increasingly apparent with subsequent proofs and generalizations of the 
Lee- Yang theorem |3| [5J [71 [IT] . Indeed, all these boil down to questions regarding 
linear operators preserving non-vanishing properties of multivariate polynomials 
and entire functions. In an attempt to provide a common ground for Lee- Yang 
type problems in statistical mechanics and various contexts in complex analysis, 
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probability theory, combinatorics, and matrix theory, we consider here the question 
if one can characterize such operators. 

Given an integer n > 1 and £1 C C" we say that / S C[z±, . . . , z n ] is Vl-stable 
if f(zi, . . . , z n ) ^ whenever (zi,...,z n ) E f2. A K-linear operator T : V — > 
K[zi, . . . , z n ], where K — M or C and V is a subspace of K[zi, . . . , z n ], is said to 
preserve fi- stability if for any f2-stable polynomial / E V the polynomial T"(/) is 
either ^-stable or T(f) = 0. For k = (ki, . . . , n n ) £ N™ let K K [zi, . . . , z n ] = {/ £ 
K[zi, . . . , z„] : deg^. (/) < k^, 1 < z < n}, where deg z . (/) is the degree of / in Zi- If 
f cC and = ^ n then O-stable polynomials are also referred to as ^-stable. 

Problem 1. Characterize all linear operators T : K K [zi, . . . , z n ] — > K[zi, . . . , z n ] 
that preserve £l-stability for a given set C C n anrf k £ N™ . 

Problem 2. Characterize all linear operators T : K[zi,...,z„] — > K[zi,...,z„] 
i/iai preserve fl-stability, where il are prescribed subsets ofC n . 

For ?i = 1, K = E, and = {z £ C : Im(z) > 0} Problems Q] [2] amount to 
classifying linear operators that preserve the set of real polynomials with all real 
zeros. This question has a rich history going back to Polya-Schur [9] and remained 
unsolved until very recently. We answer Problems Q] [2] when fi = f2i X • • • X O n and 
the f2i's are open circular domains in C, that is, open discs, exteriors of discs, or 
half-planes. We also state a general composition theorem and Grace type theorems 
for multivariate polynomials. Applications and full proofs may be found in [2]. 

2. Classification of Linear Operators Preserving Stability 

Let {Ci}™ =1 be a family of circular domains and k = (ki, . . . , n n ) £ N™. It is 
natural to consider the set A/" K (Ci, . . . , C„) of C\ x • • • x C„-stable polynomials in 
C K [zi, . . . , z n ] that have degree Kj in Zj whenever Cj is non-convex. Let also 

AA(d,...,C„) = |J AA K (Ci,...,C„). 

Lemma 2.1. Suppose that C\, . . . , C n , D\, . . . , D n are open circular domains and 
k = [k\, . . . , K n ) £ N". Then there are Mobius transformations 

C ^ 4>i{C) = Q% » + j' ; - fri c i = 1, 1 < i < n, (2.1) 

smc/i i/iai t/ie ^ invertible ) linear transformation $ K : C K [zi , . . . , z n ] — > C K [21 , . . . , z„] 
defined by 

$ K (/)(Z 1 ,...,Z„) = (Ci^+di)" 1 •■•(c„z„+d n ) K "/(0i(^i),--->n(^)) (2.2) 
restricts to a bijection between Af K (Ci, . . . , C n ) and N K (D\, . . . , D n ). 

In what follows the open unit disk is denoted by D and open half-planes bordering 
on the origin by H e = {C £ C : Im(e w () > 0} for some 9 £ [0,2vr), so that H is 
the open upper half-plane while is the open right half-plane. 

The following theorem answers a more precise version of Problem [1] for K = C 
and f2 = C\ x • • • x C n , where the CVs are arbitrary open circular domains in C. 
An analogous result (Theorem 1.2 in [U I]) solves Problem Q] for K — R, fl = Hq. 



LEE- YANG PROBLEMS 



3 



Theorem 2.2. Let n = («i, . . . , K n ) G N n , T : C K [zi, . . . , z n ] — > C[zi, . . . , z„] &e 
a linear operator, and Ci — (^~ 1 (Ho), where fa, 1 < i < n, are Mobius transfor- 
mations as in (|2.ip suc/i i/iaf £/ie corresponding <I> K defined in (|2.2p restricts to a 
bisection between J\f K (M.Q, . . . ,Ho) and N K {C\, . . . , C„) fc/. Lemma [K7]) . Then 

T : 7V«(Ci, ...,C B )-» AA(Ci, . . . , C n ) U {0} 

«/ and onZj/ z/ either 

(a) T /ias range of dimension at most one and is of the form 

T(f) = a(/)P, 

where a is a linear functional on C K [z\, . . . , z n ] and P is a C\ x • • ■ x C„- 
stable polynomial, or 

(b) i/ie polynomial in In variables z\, . . . , z n , w\, . . . , w n given by 



T 



\ \ {{a l z l + b i )(c l w l + di) + (a l w i + b l )(c i z t + di)) K " 



is Ci x • • • x C n x C\ x • • • x Cn-stable. 

The polynomial in Theorem l2.2l fb) is called the algebraic symbol ofT with respect 
to C\ x • • • x C n . 

Remark 2.1. Set z = (z\, . . . , z n ) E C n and w = (wi, . . . , w n ) £ C™. Note that if 
Ci = D, 1 < i < n, or Ci — C \ D, 1 < i < n, the algebraic symbol of T becomes 
a constant multiple of T[(l + zw) K ], while if d = Mg, 1 < i < n, 9 G [0, 27r), it 
is just a constant multiple of T[(z + w) K ], where (z + w) K = nj=i( z j + w j) Ki an( ^ 
(l + zw) K = rij=i (•"■"!" z j w j) Kj ■ If = § it is often more convenient (but equivalent) 
to choose the symbol T[(l + zw) K ]. 

Given a linear operator T : C[z±, . . . , z n ] — > C[z\, . . . , z„] we define its transcen- 
dental symbol with respect to Hg to be the formal power series 



where a! = a\\ ■ ■ ■ a n \ and z ■ w = z\W\ + . . . + z n w n with z = (z\, . . . , z n ) and 
w = (wi, . . . ,w n ). The solution to Problem [2] for K = C and fl — H l is as follows. 

Theorem 2.3. Let T : C[z%, . . . , z n ] — > C[zi, . . . , z n ] 6e a linear operator. Then T 
preserves Ho- stability if and only if either 

(a) T has range of dimension at most one and is of the form T(f) = a(f)P, 
where a is a linear form on C[z±, ... , z n ] and P is a Mo-stable polynomial, 
or 

(b) T[e~ z ' w ] is an entire function which is the limit, uniformly on compact sets, 
of 'Ho- stable polynomials. 

The analog of this result given in Theorem 1.4 of [5J I] answers Problem [2] for 

K = M and ft = HjJ. 

Remark 2.2. Note that by suitable affine transformations of the variables one gets 
from Theorem [231 a solution to Problem [5] for K = C and ft = Hi X • • • X H n , where 
the Hi's are any open half-planes in C. For instance, the analog of Theorem 
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(b) for the open right half-plane H.j is that the transcendental symbol of T with 
respect to ffi, i.e., the formal power series 



T[e z - W ] := J2 nn^J e C[zi, ...,z n ] \[w u 



w 

,w n \ 



a< 

is an entire function which is the limit, uniformly on compact sets, of -stable 
polynomials. 

3. Linear Operators in Lee- Yang Type Problems 

As explained in Part B of [2 II] , a common feature of the key steps in existing 
proofs and generalizations of the Lee- Yang theorem (Theorem 1 1 . 1 [) and Heilmann- 
Lieb theorem [3] is that they use certain linear operators on multivariate poly- 
nomials and one has to show that these operators preserve appropriate stability 
properties. Therefore, Theorems 12.21 |2~31 provide in particular a simple way of de- 
riving these key steps. We will now illustrate this with a few examples starting with 
a short proof of Theorem 1 1 . 1 1 based on the ideas in {7J combined with Theorem 12. 31 

Proof of Theorem \l.l\ Note that (b) follows from (a) by symmetry in a i— > —a. To 
prove (a) define M to be the set of functions fx : { — 1,1}" — > C whose Laplace 
transform 

<TG{-1,1}" 

is the limit, uniformly on compact sets, of H*l -stable polynomials. 

Claim: Let i, j € [n] :— {1,. .. ,n} and Jij > 0. If /i G A4 then fl^j G A4, where 
- , s \e J ^fi(a) if o~i = o-j, 

ye J,3 ^((T) if 0~i f= CTj. 

Let us show that the claim implies the theorem. Indeed, if /io : {— 1, 1}™ — > C 
is such that n(a) = 1 for all a G { — 1, 1}" then its Laplace transform Z Mo equals 
(e hl + e~ hl ) ■ ■ ■ (e' 1 " + e~ h "). Since 2cosh(x) = lim fe ^ 00 [(l + x/k) k + (1 - x/k) k ] 
and |(1 — C)/(l + C)| < 1 whenever Re(C) > it follows that fi G M. Then by 
successively applying to the transformations defined above for all pairs G 
[n] x [n] one gets (a). 

To prove the claim note that Zf lij = T(Z fla ), where 

d 2 

T = cosh(Jy) + smh(Jy)-— — . 

O Z% (J Zj 

By Theorem 12.31 and Remark 12.21 the operator T preserves Hi -stability. Since T is 
a second order (linear) differential operator, by standard results in complex analysis 
we have that if fk — > / uniformly on compacts then T(fk) — ► T(f) uniformly on 
compacts. This proves the claim. □ 

Further illustrations of the aforementioned philosophy are as follows: 
(i) Many known proofs of the Lee- Yang theorem are based on Asano contrac- 
tions or variations thereof [TJ [TT] . Let 

f(zx, . . . , z n ) = a(z 3 , ...,z n ) + b(z 3 , . . . , z n )zi + c(z 3 , z n )z 2 + d(z 3 , ■ • • , z n )z 1 z 2 

be a polynomial in n > 2 variables with deg z .(/) < 1 for i = 1,2. The 
Asano contraction of / is A(/)(zi, . . . , z n ) = a{z 3 , . . . , z n )+d(z 3 , . . . , z n )z\. 
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The key fact used in the aforementioned proofs is the following property of 
Asano contractions. Let k = (/si, . . . , K n ) G N n with n > 2 and Ki = K2 = 1. 
Then the linear operator A : C K [zi, . . . , z„] — > C K [zi, . . . , z„] preserves ID- 
stability. This is an immediate consequence of Theorem l2.2l and Remark |2. II 
since the algebraic symbol of A with respect to ID™, i.e., 

A[(l + zti)) fi ] = (1 + zw)^ 3 '-'^ (1 + ziwiw 2 ) 

is a ID-stable polynomial, 
(ii) In [7] Lieb and Sokal generalized Newman's strong Lee- Yang theorem [§]. 
A key ingredient in Lieb-Sokal's proof is the following result which they 
obtained in the process. Let {Pi(u)}™ 1 and {Qi(v)}^L 1 be polynomials in 
n complex variables u — (tti, . . . , u n ), respectively v = («i, . . . , v n ), and set 

u) = £ P<(«)Qi(«), 5(2) = Pi(d/dz)Qi(z), 

j=l i=l 

where z = (zi, . . . ,Z n ), d/dz = (d/dzi, . . . ,d/dz n ). If R is H?' l -stable then 
S is either IE -stable or identically zero. A simple proof is as follows. Define 
a linear operator T : C[ui, . . . , u n , vi, . . . , v n ] — > C[ui, . . . , u n , «i, . . . , v n ] by 

and extending linearly. The above statement is equivalent to proving that 
T preserves H.j -stability. By Theorem 12.31 and Remark 12.21 this amounts 
to showing that the transcendental symbol of T with respect to HE , i.e., 



a,f3 i=l 



is an entire function which is the limit, uniformly on compact sets, of Un- 
stable polynomials. This is clearly true since e xy = lim / ! C ^ 00 (l + xy/k) k . 
(iii) In 5 J Hinkkanen established the following composition theorem. Let f,g G 
C ( in)[zi,...,2 n ], where (1™) = (1, . . . , 1) G N n . If f,g are ©-stable then so 
is their Hadamard-Schur product (or convolution) 

(/ • 9)(z) = / (a) (0).9 (a) (0)z Q , z = ( Zl , z„), 



unless / • g = 0. This follows readily from Theorem 12.21 and Remark 12. II 
Indeed, fix a D-stable polynomial g G Cnn\[z%, . . . , z n ] and let T be the 
linear operator on CnnJzi, ... ,z n ] given by T(f) — f • g. The algebraic 
symbol of T with respect to ID™, i.e., 



T 



(1+ZW7) (1) = g(ziWi, . . . , Z n W n ) 



is obviously ID-stable, which yields the desired result. In [5] this was then 
used to argue as follows. Let et^ G C with \dij\ < 1, 1 < i, j < n, and set 

fij(zi, ... ,z n ) = (1 + OijZi + Wfz.j + ZiZj) JJ (1 + Zk), 1 < * <3 < n. 
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It is not hard to see that fij is ID-stable and by taking the Hadamard-Schur 
product of all these polynomials one gets 

(/^•••••/(n-i^o) = e zS nn a J 

SC{l,...,n} i£S j£S 

which is D-stable by the above composition theorem. This proves a strong 
version of what is usually referred to as the Lee- Yang "circle theorem" . 
(iv) Let us finally consider the proof of the Heilmann-Lieb theorem given in [3] . 
Given a graph G — (V,E), \V\ = n, equipped with vertex weights {zi}n=v 
and non-negative edge weights {A e } e e-E form the polynomial 

F G {z,X) = Yl (1+XeZiZj). 

e={i,j}eE 

Define also the linear operator MAP : C[z%, . . . , z n ] — > C^n) [zx, . . . , z n ] that 
extracts the multi-affine part of a polynomial, i.e., if f(z) = X^aeN™ a ( a ) za 
then 

MAP (/)(*) = E a{a)z a . 

a: ai <1, l<i<n 

Clearly, Fq{z, A) is H.j-stable in the z^s. Since the multivariate Heilmann- 
Lieb polynomial is given by MAP[Fg(z, A)], in order to prove the Heilmann- 
Lieb theorem it is enough to show that MAP preserves Hi -stability. Now 
the transcendental symbol of MAP with respect to Wi (cf. Remark ing is 

a n 

E ^ = IK 1 + ^)- 

a: ai<l. l<i<n i— 1 

This polynomial is -stable and thus Theorem l2.3l and Remark 1 2 . 21 imply 
that MAP preserves -stability, as required. 

4. Multivariate Master Composition Theorems and Apolarity 

The next theorem extends to several variables the Hadamard-Schur convolu- 
tion results of Schur-Malo-Szego, Walsh, Cohn-Egervary-Szego, de Bruijn [TO] and 
provides a unifying framework for multivariate generalizations of all these results. 

Theorem 4.1. Let k 6 N n and f, g G C[zi, . . . , z n , Wi, . . . , w n ) be of the form 
f{z,w) = E (f)^M^, g(z,w) = E (^)Qa(z)w a , 

a<K ^ ' a<K ^ ' 

where z = (zi, z n ), w = (wi, . . . , w n ). 

(a) If f and g are Mg -stable for some < 9 < 2n, then the polynomial 

£ (:)p.H«„H 4 £ g(o,») ■ £S(=,o) 

is Mg -stable (in 2n variables) or identically zero. 

(b) // / and g are Bi-stable, then the polynomial 

is Instable (in 2n variables) or identically zero. 
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Proof. Suppose that /, g are as in part (a) of the theorem. Let 

T : Cp[zi, . . . , z n ] -> C K [zi, . . . , z n \ and S : C K [z ± , z n ] -> C 7 [zi, . . . , z n ] 

be the linear operators whose algebraic symbols with respect to MS (cf. Remark 2.1) 
are /, respectively g, with /?, 7 £ N™ appropriately chosen. By Theorem 1 2 . 21 both 5 
and T preserve He-stability, hence so does their (operator) composition ST whose 
symbol is precisely the polynomial in (a). Applying Theorem 12 . 21 again we conclude 
that this polynomial is He-stable unless it is of the form A(z)B(w) for some poly- 
nomials A and B. If this is the case and these polynomials are not identically zero 
then A(z) must be Hg-stable (being the polynomial P in Theorem 12.21 (a)) and by 
exchanging the roles of / and g we get that B(w), thus also A(z)B(w), must be 
Hg-stable. This proves (a). Part (b) follows similarly. □ 

Remark 4.1. A still more general composition theorem (involving polynomials in 
An variables) is given in Theorem 3.3 of [21 II]. 

For two polynomials f,g£ C[zi, . . . , z n ] and n £ N™ define 

{/,<?}* ■■= £(-i)V (a) (o)<7 (K - a) (o) 

and call / and g apolar if they both have degree at most k and {/, g} K = 0. The 
classical Grace apolarity theorem [10] may be stated as follows. Let n = 1, C be a 
circular domain in C, and /, g be univariate complex polynomials of degree k > 1. 
If / is C-stable and g is C \ C-stable then {/, g} K ^ 0. Several authors asked if 
Grace's apolarity theorem could be extended to multivariate polynomials (see, e.g., 
[5]) but the precise form of such extensions remained unclear. We have the following 
multivariate apolarity theorems for discs and exteriors of discs (Theorem I4.2j) and 
for half-planes (Theorem 14.3ft . 

Theorem 4.2. Let Ci, 1 < i < n, be open discs or exteriors of discs and f,g € 
C K [zi, . . . , z n ], where n — [k%, . . . , K n ) £ N n . Suppose that 

(i) / is C\ x • • • x C n -stable and deg z .(f) — Kj if Cj is the exterior of a disk, 
and 

(ii) g is (C\ C\) x • • • x (C \ C n )-stable and deg z .(/) = Kj if Cj is a disk. 
Then{f l9 } K ^0. 

Theorem 4.3. Let C\ and C2 be two open half-planes with non-empty intersection, 
k £ N™, and f,g G C K [zi, . . . ,z n ]. If f is C\-stable, g is C^-stable, and K < a + [3 
for some a £ supp(/),/3 € supp(g), then {f,g} K ^ 0. 

Theorems 14 . 2l4"4~3l provide in particular an answer to a question of Hinkkanen [5]. 
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